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The phenomenon of random intensity patterns, for waves propagating in the presence of disorder, 
is well known in optics and in mesoscopic physics. We study this phenomenon for cold atomic 
gases expanding, by a diffusion process, in a weak random potential. We show that the density- 
density correlation function of the expanding gas is strongly affected by disorder and we estimate 
the typical size of a speckle spot, i.e., a region of enhanced or depleted density. Both a Fermi gas 
and a Bose-Einstein condensate (in a mean-field approach) are considered. 
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I. INTRODUCTION 

A wave propagating in a random medium undergoes 
multiple scattering from inhomogeneities. As a result, a 
complicated, highly irregular, wave intensity pattern is 
formed. It is characterized in statistical terms, with the 
help of various correlation functions and probability dis- 
tributions. Such intensity patterns have been thoroughly 
studied for electromagnetic waves in disordered dielectric 
media (optical speckles) [l| . They are also well known for 
the Schrodinger waves in disordered electronic conduc- 
tors, where, due to the electron scattering on impurities, 
a random distribution of electron density and currents 
is established within the sample. These "Schrodinger 
speckles" manifest themselves in various mesoscopic ef- 
fects, for instance, in the sample-to-sample conductance 
fluctuations 

Cold atomic gases offer a new playground for studying 
matter waves in the presence of disorder. Unlike the case 
of electrons in disordered conductors, for atomic gases 
one can directly measure the atomic spatial density or, 
more precisely, density integrated along the direction of 
imaging. There have been a number of recent experi- 
ments S2i,i3i, as well as some theoretical work 
[1, Q on the effect of disorder on atomic Bose-Einstein 
condensates (BEC), in quasi-one-dimensional geometry. 
Some of these experiments reported observation of appar- 
ently random albeit reproducible interference patterns - 
the atomic analog of the Schrodinger wave speckles in 
disordered conductors. 

In the present paper, we develop a theory for such 
"atomic speckles" in two and three spatial dimensions 
when, unlike the one-dimensional (ID) case, atoms can 
propagate by diffusion. The diffusion of the BEC in a 
random potential has been recently studied p^ . [Til . [T^ 
and it was pointed out that only some fraction of the con- 
densate diffuses away, while the rest stays localized near 
its initial location (Anderson transition for BEC) (Tll.[l^. 
Only the average density of the evolving atomic cloud was 
studied in Refs. [O, The atomic speckles, however. 



are characterized by more complicated quantities, such 
as the density-density correlation function studied in the 
present paper. Before plunging into calculations, let us 
emphasize that we assume no (or negligibly small) dis- 
order within the trap, i.e., the atoms are released from 
a "clean" trap and diffuse in the random potential. (A 
different setup is when the disorder, created within the 
trap, has a strong effect already at the equilibrium state 
of the system, causing a transition to an insulating phase 
p^.) We focus mainly on the case of Fermi gases (Sec- 
tion |Tl|. A brief discussion of the BEC case is given in 
Section Hm 



II. 



ATOMIC SPECKLES IN A DIFFUSING 
FERMI GAS 



We consider N fermions, initially occupying the N low- 
est energy eigenstates in a trap. Assuming zero temper- 
ature and neglecting interactions we can write the quan- 
tum expectation value of the particle density as 



(n(r))o = ff.V (r)r 



(1) 



where — 25' -1-1 is the spin-degeneracy factor and 0„(r) 
is the orbital part of an eigenstate in the trap. The factor 
/„ is unity for an occupied site and zero otherwise. 

At time t = the Fermi gas is released from the trap 
and the single-particle wave functions start to evolve ac- 
cording to the Schrodinger equation 



ihdMv,t) = -— AV„(r,i) + V^(r)^Jr,i), (2) 



where V{v) is the external static random potential. The 
quantum expectation value of the particle density, at 
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time t, for a given realization of randomness is equal to 

n 

= g, jd^R Jd'^R' Gir,-R,t)G*ir,-R',t) 



(3) 



where G(r, R, t) is the retarded Green's function of the 
Hamiltonian in Eq.(l2])- The sum in Eq.Q can conve- 
niently be written as 

Ep 

= -\ j deImg^(R.',R,£), (4) 

— oo 

where e„ are the particle eigenenergies in the trap and 
is the Fermi energy. The Green's function tJ(R', R, e) 
refers to a particle in the trap and it should not be 
confused with the function G{y, R, t) which describes 
propagation from point R to an observation point r, 
upon the release of the gas from the trap. 

The density (ri(r, t)) is still a random variable, in the 
sense that its value depends on the specific realization 
of the random potential V{v). Denoting averaging over 
realizations by an overbar, we write 

(n(r,t)) = jd'^R j d'^R' G{v,K,t) G* (r,R',t) 

Ep 

X J delmg{R\R,e). (5) 

— oo 

The last integral in this expression contains information 
about the initial state of the gas in the trap. The product 
of the two propagators, GG*, propagates this informa- 
tion in space and time. 

In experiment, an ensemble of many random samples 
can be prepared either by creating different realizations 
of the random potential or by making slight changes in 
the Fermi energy ep. It is well known in mesoscopic 
physics that, for a fixed realization of randomness, vari- 
ous properties of a system are extremely sensitive to the 
precise position of Sp Indeed, an "ensemble of disor- 
dered conductors" is often produced from a single sample, 
whose randomness is fixed by technology, by changing 
the Fermi energy. In this respect, cold atomic gases con- 
stitute a more convenient system for studying disorder 
related effects. A different realization of disorder can be 
easily created by changing the optical speckle pattern, 
on which the atoms are being scattered. Furthermore, 
since the precise value of £p cannot be controlled, from 



one experimental run to another, to an accuracy higher 
than a few percent it appears that in experimenting with 
cold atoms an "ensemble of random samples" is created 
automatically, in a natural way. 

The average product of the two Green's functions 
in Eq.([ni) is a standard object in the theory of disor- 
dered systems. In the diffusion regime, considered in the 
present paper, it can be transformed as [l^ 



G(r, R,t) G* (r, R',i) 
1 



dePjr,R,t)ImG(R-R',e), (6) 



where 



P,(r,R,t) 



,-|r-R|V4D,t 



(7) 



is the diffusion propagator, in d dimensions (d = 2, 3) 
and is the diffusion coefhcient for a particle at energy 
e. The average Green's function is given by [ij 

G(Ai?,e) = Go(Ai?,£)e-^-^/2'= , Ai? = |R-R'|,(8) 

where Gg is the free Green's function and is the par- 
ticle mean free path. The Eqs. ([7|) and ^ for P and G 
are valid only for sufficiently large values of the energy 
parameter e, namely if A:^?^ ^ 1, where fc^ = a/ 2me / It? . 
Therefore, kplp is the essential parameter which deter- 
mines the overall behavior of the gas, after switching off 
the trap, where fcp and Ip denote the value of the wave 
number and of the mean free path at the Fermi energy 
ep. Only for fcpZp 3> 1, which is the case assumed in 
this paper, most of the atomic cloud will diffuse away 
from the trap. (A similar condition, with the chemical 
potential fi replacing the Fermi energy, is required for the 
diffusive behavior of a BEG cloud 

To facilitate analytic treatment, we assume that the 
size of the trap is much smaller than its distance from 
the observation point r. Then, choosing the coordinate 
origin somewhere inside the trap, one can set R = in 
the second argument of the diffusion propagator, writing 
it simply as P^{r,t) [r ~ |r|). Moreover, since upon 
release from the trap, each particle goes on its own (long) 
diffusive trajectory, the actual shape of the trap is of 
no importance. It is convenient to replace the actual 
harmonic trap by a cubic trap of size L, with periodic 
boundary conditions. Finally, Fourier transforming G in 
([5]) with respect to R — R', and performing integration 
over R and R' in ([5]), we obtain 

= ldePAr,t)ImG{Ke) 



Cp 

X J de'lmg{k,e') . 



(9) 



Since we assume no disorder in the trap, it follows that 
Imt/(k, e') = —TrS{e' — e^) and integration over e' results 
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in a step function 0(fcp — k). Furthermore, since fcpZp ^ 
1, it follows that the weak disorder condition, k^l^ ^ 
1, is satisfied for the great majority of fc's in Eq. ([9]). 
Therefore, ImG(k, e) can be approximated by — 7r(5(e— e^.) 
and integration over e can be carried out, resulting in 



{h{r,t)) = g.L" 



(10) 



where the diffusion kernel Pf.{r,t) is given by Eq. ([7]), 
with R = and e set equal to h^k^ /2m. Equation 
(llOp has a very simple interpretation: Particles, prior 
to their release from the trap, occupy all states up to 
fcp, with L'^d'^k/{2TT)'^ being the number of particles 
in the element d'^k. When the trap is switched off, 
particles start diffusing and P^, (r, t) is the probability 
density that a particle with wave number k will reach 
point r in time t. Integration over k gives the average 
particle density t)). An equation similar to (jlOp 
exists also for a BEC It was emphasized there, and 
the same is true for the present case of a Fermi gas as 
well, that some fraction of the released particles will not 
propagate by diffusion but will remain localized near 
their original location. For kplp ^ 1 the fraction of such 
particles is small and we ignore them in the present work. 

We turn now to the calculation of the density-density 
correlation function. Let us first note that, already in 
the absence of disorder, a Fermi gas possesses some sub- 
tle density correlations of purely quantum nature p^ . 
The expectation value, (n(r,t)), is ignorant about these 
correlations. However, as has been particularly empha- 
sized in [l^, [l^ , a single imaging experiment (with suf- 
ficient resolution) does not measure (n(r,t)) but rather 
one particular event, i.e., some particular density pat- 
tern, n(r,t), whose probability is dictated by the many- 
body wave function of the system. The density correla- 
tion function can be extracted from such noisy density 
patterns by averaging the product n(r)n(r') over many 
experimental realizations, taken under (as far as possi- 
ble) identical experimental conditions, for fixed r, r' and 
t. Such averaging would yield the theoretically calcu- 
lated quantum expectation value (n(r)n(r')). It should 
be noted, however, that already a single experimental im- 
age, although noisy and "grainy", contains information 
about the density correlation function. This information 
will be revealed by taking the product n(r)ri(r + Ar) and 
averaging it over many points r, for fixed Ar (the equiv- 
alence of the two averaging procedures constitutes the 
"ergodic assumption" in mesoscopic physics [J). 

For fermions confined to a trap of size L, the density- 
density correlation function exhibits rapidly decaying os- 
cillations with a characteristic period Axq « fcp^ ~ 
LN^^/'^ [i3,[i3- When the gas is released from the trap 
it starts expanding, and, for times t > (L/v-p) = tg, the 
size of the cloud grows linearly with time and so does the 
correlation length Ax{t) w Aa;Q -t/tQ [ll|- Thus, roughly 
speaking, the free, ballistic expansion amplifies the scale 



of correlations by the factor t/tQ. 

Below we show that in the presence of a random po- 
tential, i.e., when the expansion is diffusive instead of 
ballistic, the picture is different: The size of the atomic 
cloud grows as Vt whereas the short-range correlations 
do not get amplified at all. The density-density correla- 
tion function in the presence of disorder is defined as 



C(r,r',t) = (n(r,t)n(r',i)) - {n{r,t)) {h{r' ,t)) 



-(5(r-r')(n(r, t)) 



(11) 



The last term describes trivial correlations, which exist 
already in a classical ideal gas and which are commonly 
subtracted, in order to isolate the nontrivial correlations 
[3] . There are two kinds of averaging in (fTT|) : The quan- 
tum mechanical averaging, for a given realization of dis- 
order, and averaging over the ensemble of different real- 
izations. The first averaging is straightforward and leads 
to [14]: 

(n(r,t)n(r',<)) = 52^/„K(r,t)|2^/„J^„(r',i)|' 



^/„C(r,t)V„(r',t) 



+ gJ{r-r')J2fn\i^nir,t)\\ (12) 

n 

where the absence of overbars indicates that this expres- 
sion refers to a specific realization of the random poten- 
tial. Next, we must average ([T2|) over the disorder. This 
involves averaging products of four single-particle wave 
functions. For short-range correlations, i.e., on a scale 
smaller than the mean free path , such averages decou- 
ple into products of pairwise averages [ij. For instance. 



7/;*(r,^)^„(r^^)fc(r^t)7^„(r,t) 
« ^;(r,i)V'„(r',t) <(r',i)V'™(r,t) 



+ ^;(r,i)V'™(r,t) <(r',i)V„(r',t) 
Performing such decoupling in Eq. (I12p and subtracting 



r-r'l < /. 

(13) 



the product of averaged densities, (n(r,t)) (n(r',i)), we 
obtain 

C{r,r',t) = ^/„/„ |-g^A„„(r,r',t)A,*„„(r,r',t) 



~ 9s ^nm i)^nm (l"' i 1"' J *) 

+ g2|A„,„(r,r',i)|'} , 



where 



A„„,(r,r',t) EE V;(r,t)V'™(r',i), 



(14) 



(15) 



The first term in p4p describes quantum correlations, 
due to the Pauli exclusion principle. In particular, for 

2 

r' r, it approaches the value —g {h{r,t)) and, thus. 
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it is proportional to 7V^. The third term is of "classical" 
origin, in the sense that it originates from the interference 
between multiply scattered waves. It contributes positive 
correlations, similarly to speckle pattern in optics. How- 
ever, in contrast to a single frequency laser speckle, here 
there are many waves with different frequencies. Since 
contributions from different frequencies should be added 
incoherently, i.e., intensities (rather than amplitudes) are 
summed up, the third term is proportional to N and it 
will be neglected. The second term in ([H]) is a com- 
bination of quantum and classical correlations. Its sign 
and the factor originate from the exclusion principle. 
Since, however, wave functions at different energies are 
essentially uncorrelated, this term is also proportional 
to N and can be neglected in the large- iV limit. Thus, 
keeping only the first term in (|14p . we obtain 



For d = 3 [H, 



C(r,r',t) = -g. 



-.gjF(r,r',i)r 



(16) 



The function F(r, r',i) = /nV'n(r, i)V'n(r', i) can be 

n 

expressed in terms of the Green's functions, in complete 
analogy with the earlier derivation for the average den- 
sity. The only difference is that now there are two "ob- 
servation points", r and r'. The resulting expression for 
F(r, r',i) (compare to Eq. (O) is 

F(r,r',t) = Jd'^R J d'^R' G{r, R,t) G* {r',R\t) 



X J delmg{R\R,e) 



(17) 



In order for the average product of Green's functions not 
to be exponentially small, the pair of points r, r' should 
not be separated by more than the mean free path. On 
the other hand, |r — R| is much larger than /p. The com- 
putation of the average product in p7p is straightforward 
and the result is 

G(r,R, t) a* (r',R',t) = -i J de f{Ar,e) P, {r,R,t) 

xImG(r-R',e), (18) 

where the extra factor f{Ar, e), as compared to Eq. ©, 
is given by 

f(Ar,e) = —lmG(Ar,e), Ar = |r-r'| . (19) 

Following the same line of derivation as for the average 
density, Eq. (fTUl). we arrive at 

F{ry,t) = I ^P,{r,t)fiAr,e,). (20) 



sin(fcAr) _Ar/2/, 
kAr 



(21) 



While the kernel Pk and the mean free path Ik are slowly 
changing with fc, the function f{Ar,ek), for Ar kp^, 
contains a rapidly oscillating factor within the integration 
region in (PO)) . Taking the "slow" functions out of the 
integral and computing the remaining integral we obtain 

C(r,r',t) = -g^L''^P^^ir,tre-^^/^^^ 

^ Jsin(fcpAr) ~ (fcpAr) cos(fcpAr)]^ 
(fcpAr)*^ 

where Pj. and Zj, are taken at fc = fc^. It is convenient to 
normalize the correlation function by the average density 
which yields 



/(r,r',t) 



C{v,v',t) 



{h{r, t)){n{r',t)) 



9s^ 



9 -Ar/i„ [sin(fcFAr) - (fcpAr) cos(fcFAr)]' 



(fcpAr) 



(23) 



The analogous calculation for d = 2 yields 



/(r,r',0 = -le--/'^^il^, (24) 
9s (fcpAr) 



where J]^(fcpAr) is the Bessel function of the first kind. 

The normalized density-density correlation functions, 
given m Eqs. (HSl) and JM]) are plotted in Fig. [1] (dashed 
lines). The solid lines in this figure represent the cor- 
responding functions obtained by computing the inte- 
gral in Eq. numerically, without taking the "slow" 
functions out of the integral. In this computation a 
white noise random potential has been used which yields 
Zj. = const in three dimensions and If^ ^ k in two dimen- 
sions (this implies that the diffusion coefficient is propor- 
tional to k and fc^ in three or two dimensions, respec- 
tively). Figure [T] demonstrates that for the decaying en- 
velope of /(Ar) the agreement between the numerically 
exact results and the approximate expressions, Eqs. (j23|) 
and ([M)) . is quite good. For the oscillations the agree- 
ment is only qualitative. Due to the rapid decay of the 
envelope function, the oscillations are rather small and 
clearly visible only on an amplified scale, as shown in the 
insets to the figure. 

It is quite remarkable that, while the atomic cloud 
keeps expanding, the local normalized density correla- 
tions, Eqs. ((23|) and (|24|) . do not depend on time. 
In particular, the characteristic length of oscillations, 
Ar ~ fcp^, remains the same as for the gas in the trap, 
prior to its release. This behavior is in sharp contrast 
with that for a free ballistic expansion, when the spa- 
tial oscillation period is growing linearly in time, just as 
the average interparticle distance. The essential differ- 
ence between ballistic and diffusive cases can be traced 
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FIG. 1: The normalized density- density correlation function 
I{Ar,t) of a Fermi gas with Qs — 2, for kplp = 10, r — 50Zf 
and t — 3r^ /{2dDkp)- (a) corresponds to d = 3 and (b) to 
d = 2. The dashed curves in (a) and (b) correspond to Eqs. 
(|23|) and (|24[) . respectively. The solid curves are the result of 
an exact numerical integration, as described in the text. The 
insets show the oscillating decay of I{Ar,t) with increasing 
Ar. 



to the evolution of the phase of a wave function. In the 
ballistic expansion the characteristic spatial period of a 
wave packet keeps increasing. (This statement can be 
easily verified by taking an initial wave packet of a Gaus- 
sian shape, 0(r) ~ exp(— ^ -I- ip-), with b <ti a, whose 
time evolution is exactly soluble.) On the contrary, in 
the diffusive expansion the phase of a wave function gets 
randomized after few scattering events and no long-range 
order in the phase can be established. An oscillating wave 
function with some period b, subjected to diffusive evo- 
lution, will locally look as a plane wave, with the same 
period, where "locally" means on a scale smaller than 
the mean free path. Thus, at any time t, the diffusing 
wave function can be viewed as made up of " patches" of 
plane waves, of size each, but with no phase relation 
among different patches. This observation explains the 
somewhat counterintuitive behavior of the correlations, 
namely, that the correlation functions in Eqs. and 
(j24p remain stationary, while the interparticle distance 
increases under the expansion. This leads to an increase 



of the relative fluctuation of the particle number, in a 
given volume. For a homogeneous Fermi gas in equilib- 
rium (in three dimensions) the particle number variance 
AiV^, in a certain volume, is not equal to the average 
number of particles N (in the same volume) but is pro- 
portional to iV^ ^ log [T^H^]- This means that, due to 
correlations, the Fermi gas possesses some kind of "rigid- 
ity". (The effect is particularly spectacular in ID, where 
AA^^ grows only as log A). Free expansion of the gas, 
in the absence of disorder, does not affect this rigidity, 
because, as was already mentioned, the scale of correla- 
tions is amplified in exact proportion to the interparticle 
distance. The disorder disrupts this proportionality and 
leads to destruction of rigidity and to the AA^ — A be- 
havior (in the long time limit). Therefore, the image of a 
Fermi gas, expanding in the presence of disorder, should 
look more " grainy" than the image of a freely expanding 
gas. This effect might be observable experimentally. 



III. CORRELATIONS IN A DIFFUSING BEC: A 
MEAN-FIELD APPROACH 

So far we considered the dynamics of a degenerate 
Fermi gas. We now briefly discuss the case of a BEC, 
expanding in the presence of a random potential. Within 
the mean-field approach the BEC is described by a 
macroscopic wave function, \['(r,t), whose dynamics sat- 
isfies the Gross-Pitaevskii equation 



2m 



A*(r,t) -f I/(r)'J'(r,t) 



+ g\^{r,t)\H{r,t), 



(25) 



where g is the interaction parameter related to the scat- 
tering length (we assume positive g, i.e., repulsive inter- 
actions). Equation (|25p describes the evolution of the 
condensate, in the random potential V{r), upon its re- 
lease from the trap. We assume an isotropic harmonic 
trap, characterized by frequency u. For weak random- 
ness the expansion can be separated into two distinct 
stages d, [nl, nil : a rapid ballistic " explosion" , during 
the time of order (1/w), followed by an essentially linear 
evolution. The first stage is dominated by the nonlinear- 
ity. At the second stage, however, most of the interaction 
energy had been already converted into the kinetic (flow) 
energy, so that the nonlinearity becomes weak and is ne- 
glected |2l| . Thus, below we consider the linear equation 



ihdt'^ir^t) = -—A^{r,t) + V{r)^{r,t). (26) 
2m 

The initial condition for this equation is supplied by the 
wave function <I>(r) at the end of the first stage of the 
expansion, i.e., at time of order l/w. Qualitatively, this 
wave function is of the form (23, u3| 



$(r) = T{r) exp(zr^/ao), 



(27) 



6 



where gq = (fi./mw)"'^/^ is the oscihator size of the trap 
and J^{r) is an envelope function which decays on the 
characteristic distance Rq ^ ap, where Rq is the initial 
size of the BEC in the trap. The envelope is often ap- 
proximated by a Gaussian, !F{r) = A exp{—r^ /2Rq), or 
by the inverted parabola J-{r) = Ay/1 — (r/Ro)'^, with 
A being the normalization constant. 

The emerging linear problem, Eqs. (|26|) and ((27)) . is 



considerably simpler than the fermionic problem treated 
above. This is because the BEC is described by a single 
coherent wave function which can be treated as a classi- 
cal field. Within such mean-field description, there are 
no density fluctuations in the absence of disorder. The 
notion of density correlations, with their inherent statis- 
tical features, becomes meaningful only in the presence 
of an external random potential. Thus, the problem be- 
comes similar to that considered in the theory of opti- 
cal speckles, where a classical electromagnetic wave or a 
scalar wave propagates and gets scattered on a random 
potential [Ij. The essential difference is that in the the- 
ory of optical speckles one usually assumes a monochro- 
matic field, whereas the BEC wave function, Eq. 
contains a broad spectrum of wave numbers. For a 
monochromatic field V'(^(r) with a wave number ki^, the 
disorder-induced intensity-intensity correlation function, 

C„(r,r') = |^*(r)■^/'„(r')|^ was calculated in [H, 



(28) 



a(r,r') = n^(r)n^(r') Jo'(fc^Ar)e-^'-/'-, (d=2). 



(29) 



where n^{r) — \ip^{r)\ is the average intensity of the 
wave. 

The extension to the nonmonochromatic case of a BEC 
is quite straightforward. Let us start with the average 
density of the BEC, n{r). The BEC wave function at 
time t is given by 



*(r,t) = / d'*i?G(r,R,t)$(i?). 



(30) 



Using (O, one obtains in the large t (and large r) limit 

m 



n(r) = |*(r,t)|2 

= ldePAr,t) I ^ImG(k,£)||.(fc)p(31) 

where $(fc) is the Fourier transform of the initial con- 
dition ((27|) . The (equal time) field-field correlation func- 
tion, Cfioid(Ar, t) = 5'*(r,<)*(r',i), differs from only 
by an extra factor, f{Ar,e), defined in p^ . 



(27r)° 
xImG(k,e) |$(fc)|^ 



(32) 
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FIG. 2: The normalized density-density correlation function 
/(Ar, t) of a BEC for (a) d = 3 and (b) d = 2. As in Fig. [T]the 
parameters were set to fcoZo = W, r — 50/o, t — 3r'^ / {2dDko) , 
Dk = Dkoik/kof-", Ik = lo {k/kof-". 



Approximating, as above, ImG(k, e) = —nS^e — £j,), we 
obtain the final expression for the field-field correlation 
function, 



Gfieid(Ar,t) = 



d'^k 



{2n 



^^P,ir,t)fiAr,Sk)\m\- (33) 



The short-range density-density correlation function is 
given by 



C{Ar,t) = n(r)n(r') — n(r) n(r') — (5 (r — r') n(r) 



|Gfieid(Ar,t)r 



(34) 



Equations ([33|l and p4|) , supplemented by the expression 
(fro)) for the / function, provide the general solution for 
the density correlations in a BEC diffusing in a weak 
random potential. 

To get an estimate for the three-dimensional normal- 
ized correlation function, 



/(Ar,i) = 



g(Ar,t) 
n(r', t) n(r', t) ' 



(35) 
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we assume a Gaussian envelope, ^{f) = 
A exp(— r^/2i?g). Then the Fourier transform of 
the BEG wave function, Eq. (|27p. takes the form 

||.(fc)P = AeM-k^kl), (36) 

with the normahzation constant A and fco ~ '^^o/'^o — 
1/^, where ^ is the healing length. The numerical evalu- 
ation of Eq. ([55]) (see Fig. [2]) shows that, quite similar 
to the case of the Fermi gas, /(Ar, t) decays on a scale 
Ar ^ k^^. This scale remains nearly constant in time. 
Since the "cutoff function", |$(fc)p, under the integral in 
Eq. ([M)) is not as sharp as the step function in Eq. (PO)) , 
the oscillations of the normalized correlation function are 
even weaker than for the fermionic case. The short-range 
correlations, shown in Fig. [21 imply that the image of a 
condensate, diffusing in a weak random potential, should 
exhibit a random pattern of particle density (speckle). 
The typical size of each speckle spot is of the order of a 
few healing lengths. 

IV. CONCLUSIONS 

We have studied density fluctuations in a cloud of cold 
atoms, expanding in the presence of a weak random po- 
tential. Only the simplest, two-point correlation func- 
tion, C(Ar, t), was considered. We find that the disorder 
has a strong effect on C(Ar, t), for either a Fermi gas or 
for a BEG. In both cases we obtain a random density pat- 
tern consisting of speckles of high and low density. The 
typical speckle size is determined by the decay length of 
the correlation function. For fermions this size is of the 
order of the Fermi wavelength Ap = 27r/fcp, where fcp is 
the Fermi wave number in the trap. It is interesting, and 
somewhat counterintuitive, that while the gas keeps ex- 
panding, the typical speckle size does not change. This 
is in contrast to the "clean" case of the free expansion, 
when a fixed, time-independent relation exists between 
the correlation length and the interparticle distance. For 
a BEG we confined ourselves to a mean-field treatment, 
based on the time-dependent Gross-Pitaevskii equation. 
In this case the speckle structure is caused solely by the 
random potential. The typical speckle size is given by the 
healing length of the condensate in the trap and, once 
again, this size does not change in the process of the 
expansion. If one were to go beyond the mean-field de- 
scription of a BEG, then density correlations (bunching) 
would appear even in the absence of disorder. The com- 
bined effect of quantum fluctuations and disorder on den- 
sity correlations for bosons (similarly to what has been 
done for fermions in Sec. |TT| is an interesting problem 



which is, however, beyond the scope of this paper [25|. 
The only important parameter in our theory is fcpZp, for 
fermions, and /cq^o the BEG. We have assumed zero 
temperature in all of the calculations. The extension to 
finite temperatures, at least for the Fermi case, is quite 
straightforward but has not been done in the present pa- 
per. 

Only short-range correlations were considered. It is 
well known in optics, as well as in the mesoscopic physics 
of disordered conductors, that in addition to the strong 
short-range correlations there are also weak long-range 
correlations. Such correlations manifest themselves as 
fluctuations in the transmission coefficient through a dis- 
ordered slab or as the universal conductance fluctuations 
in disordered conductors [H . In order to observe similar 
effects for cold atoms one should create a random poten- 
tial in some region (say, of a shape of a slab) and then let 
an atomic cloud impinge on that region. One could then 
take images of the transmitted, as well as of the reflected 
clouds. 

Our numerical estimates have been made for a white 
noise random potential, which is the case when the corre- 
lation radius of the potential, Rc, is smaller than the rele- 
vant wavelength. A, of the matter waves (the Fermi wave- 
length for fermions or the healing length for the BEG). 
Since the random potential for atoms is commonly pro- 
duced by creating a random pattern of light intensity (an 
optical speckle), one can assume that Rc is of the order 
of a few light wavelengths. If Rc exceeds A, one enters 
the regime of the correlated potential (colored noise). In 
that regime the diffusion coefficient is determined by the 
transport mean free path, Z^j., which is much longer than 
the scattering mean free path If^ appearing in the average 
Green's function. However, since both only weakly de- 
pend on k our results should remain qualitatively correct, 
at least in the well developed diffusion regime, i.e., when 
the cloud has spread to a distance much larger than l^^. 
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